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Abstract. We investigate the effects of competition between two complex, PT- 
symmetric potentials on the PT-symmetric phase of a "particle in a box" . These 
potentials, given by Vz(x) = iZsign(x) and V^(x) = i£ [8(x — a) — 8(x + a)], 
represent long-range and localized gain/loss regions respectively. We obtain the 
PT-symmetric phase in the (Z,£) plane, and find that for locations ±a near the 
edge of the box, the PT-symmetric phase is strengthened by additional losses to 
the loss region. We also predict that a broken PT-symmetry will be restored by 
increasing the strength £ of the localized potential. By comparing the results for 
this problem and its lattice counterpart, we show that a robust PT-symmetric 
phase in the continuum is consistent with the fragile phase on the lattice. Our 
results demonstrate that systems with multiple, PT-symmetric potentials show 
unique, unexpected properties. 



1. Introduction 

Over the past fifteen years, there has been tremendous progress in the field of parity 
and time-reversal (VT) symmetric quantum theory [TJ [2 [3] 0] ; this progress is further 
spurred by recent experiments on optical waveguides [5, 6, .7] and electrical circuits [8]. 
As it was in the case of original quantum theory, the theoretical developments have 
occurred on two distinct fronts. The first is investigations of Schrodinger equation in 
the presence of non-Hcrmitian, 'PT-symmetric, complex potentials V(x) = V*(— x) ^ 
V*(x) [3 HP] . The second is the study of non-Hermitian, PT-symmetric band 
matrices (TTJ [12l [13] that can represent the Hamiltonian for a tight-binding lattice 
with on-site, balanced, loss-and-gain potentials [HE]. 

The developments on these two fronts have been with little overlap. The 
continuum studies have focused on properties of the Schrodinger equation in the 
complex co-ordinate plane [13]; notable early exceptions are "particle in a box" 
problem with two simplest potentials, Vz{x) — iZsign(x) [TS] [THJ [T7] and V^(x) — 
it; [S(x — a) — 5(x + a)] 18]. In either of these cases, there is a nonzero threshold, 
measured in the units of relevant energy scale, for the emergence of complex 
eigenvalues. In corresponding lattice models with complex, constant potential or a pair 
of impurities, the PT-breaking threshold, measured in the units of tunneling, vanishes 
as lattice size N diverges [15] except in the case of closest or farthest impurities [2U1I2T] . 
However, on both fronts, the PT-symmetric potential is typically characterized by a 
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single parameter and the energy eigenvalues are real as long as the magnitude of that 
parameter is smaller than the threshold value. 

In this paper, we study the effect of competition between smooth (Vz) and 
localized (Vf) complex potentials on the 'PT-symmetric phase diagram. We consider 
a particle in a box of size 2L in the presence of V(x) — Vz{x) + V^(x), and its lattice 
counterpart, and calculate the phase diagram in the (Z, £) plane as a function of the 
localized-gain position a. Surprisingly, we find that a broken PT-symmetric state 
induced by one potential can be restored by additional losses to the loss region, with 
corresponding gains in the gain region, introduced by the second potential. We find 
that the phase diagrams for a < L/2 and a > L/2 are qualitatively different. With 
scaling analysis, we show that the robust phases of the continuum model are consistent 
with the fragile phases of its lattice counterpart. 



2. Particle in a box with PT-symmetric potentials 

Let us consider a particle of mass m confined to a line \x\ < L with Dirichlet boundary 
conditions for its wave function, ip(±L) — 0. When Z > 0, the complex potential 
Vz{%) = iZsiga{x) represents a constant gain potential for x > 0, with a mirror- 
symmetric loss region for x < 0. On the other hand, the complex potential V^(x) 
represents localized gain (z£) and loss (— i£) at mirror-symmetric positions ±a for 
£ > 0. The Schrodinger equation for eigenfunctions of the particle is given by 

ti 2 

- ^~d^ n (x) + V(x)ip n {x) = E n ip n (x). (1) 

It is straightforward to parametrize the eigenfunctions as [18] 

!Asin[k n (L — x)] a<x<L, 

B cos{k n x) + Ck* n s'm(k n x) < x < a, . . 

B cos(k* n x) + Ck n sm(k* n x) -a<x<0, ( ' 

Fsin[fc^(L + x)} -L<x<-a. 

This parametrization guarantees that the wave function and its first derivative are 
continuous at the origin, and that the wave function vanishes at ±L. The effect of 
(5-functions at x = ±a on the eigenfunctions is given by the discontinuity in their 
derivatives at locations ±a, 

d x M±a + ) - dM±a~) = Ti{2m{;/h 2 )ip n (±a). (3) 

Eq.([3| along with the continuity of the wave function at x — ±a provide four relations 
among the four unknown coefficients A, B, C, F. The complex wave vector k n satisfies 
E n — iZ = h 2 k 2 j2m [18]. In the following analysis, we choose ti 2 /2m = 1 = L, use the 
dimensionless quantities k n — k n L, p = a/L, Z = Z/(h 2 /2mL 2 ), and £ = t;/(h 2 /2mL) 
and employ the notation C u = cos(u) and S* = sin(w*) to simplify calculations. The 
resulting characteristic equation for k n is 

3? [\k n \ 2 k n c kn st n +e\s { i~ P )kA 2 k n c pkn s; kn } + 

«£3 |fc„| 2 S'( 1 _ p )fc n CpA; n S'fc 7i — k n S* j _ p )fc n Spk n Cfcn = 0- (4) 

We characterize the wave vector as k n = ira n + ifi n - Simultaneous solutions of Eq.Q 
and the constraint — Z = 2ira(3, obtained graphically, lead to the permitted values 
of k n and dimensionless energies E n = (na) 2 — (3 2 respectively. As the strength of 
the complex potential is increased, two adjacent or next-nearest neighbor eigenvalues 
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become degenerate and then presumably complex, leading to the breaking of the VT 
symmetry. When £ = 0, the threshold potential strength is Z c ~ 4.48 [T5] whereas 
for Z = and a = L/2, it is given by £ c ~ 5.06 [TS]. In the two cases, the eigenvalues 
E n are real for — Z c < Z < Z c and — £ c < £ < £ c respectively. When both of them 
are nonzero, the "PT-symmetric phase in the (Z, £) plane can be characterized by two 
boundary functions, Z c _(£) < Z < Z c+ (^) with £ > 0. 

It is straightforward to verify following properties of the characteristic equation: 
i) When Z = 0, the term in Eq.Q proportional to i£ vanishes, ii) When p — > 0, 1 
the ^-dependence of Eq.Q vanishes, iii) When £ = 0, the eigenvalues E n are even 
in Z. iv) The PT-phase boundaries for £ < are given by Z c+ (^) = — Z c -(— £) and 
Z c -(f) = — Z c+ {— £). Therefore, we focus only on the region £ > in the (Z,£) plane; 
when Z = = £, Dirichlct boundary conditions imply that k n — ixct n + ij3 n — mr/2. 

3. "PT-phase diagram and the re-entrant PT-symmetric phase 
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Figure 1. "PT-phase diagram for the potential Vj. The left-hand panel shows 
the threshold strength §c(p) as a function of the fractional position p = a/L of 
the localized-gain potential. Apart from a local maximum at p = 0.5 |18| . the 
threshold §c(p) diverges asymmetrically as p — > 0, 1. The center panel shows the 
quantum number n for which the energies (E n , E n +i) become degenerate at the 
VT symmetric threshold. For p < 0.5 (red triangles), the "PT-symmetry breaks at 
the ground-state energy n = 1. For p > 1/2 (blue diamonds), the corresponding 
index n(p) increases and diverges as p — > 1. The right-hand panel shows the 
scaling of the diverging threshold strength £ c as function of fractional distance 
5 <C 1 from the origin [p = 0, blue diamonds) or the boundary (p = 1, black 
circles). The vertical axis denotes l/£ c (<5) and shows that £c(<5) oc <5 -1 . 



We start this section with results for the PT-symmetric threshold £ c (p) when 
Z = 0; the properties of this model near p = 0.5 have been extensively explored |I8J. 
The left-hand panel in Fig. [I] shows the phase diagram £ c (p) as a function of < p < 1. 
Starting from the local maximum at p = 0.5, as the localized-gain location moves 
towards the origin, p — > 0, or the boundary, p 1, the threshold strength £, c (p) first 
decreases and then increases asymmetrically. The center panel shows the index n for 
which the eigenvalues E n and E n+ \ become degenerate and then presumably complex 
when the PT threshold is crossed. For p < 0.5, n — 1 (red triangles) whereas for 
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0.5 < p < 1, the eigenvalues that become degenerate have increasingly higher energy 
(blue diamonds) j|] The right-hand panel shows the scaling of the threshold strength 
for impurity locations close to the origin, 5 = p, (blue diamonds) or the boundary, 
6 = 1 — p (black circles). For 5 <C 1, we see that threshold diverges as £ c (p) oc S^ 1 
although the prefactor for the divergence is smaller when p — > (blue diamonds) than 
when p 1 (black circles). Thus, the 'PT-symmetric phase diagram at Z — is 
acutely sensitive to the distance between the localized gain and loss positions. 
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Figure 2. "PT-symmetric phase diagram in the (Z, £) plane. The left-hand panel 
shows the PT-phase boundaries Z c ±(£) for £ ^ and different locations p < 0.5 
of the localized-gain impurity. In general, as £ increases, the threshold strengths 
Z c ±, shown by dashed and solid lines respectively, decrease. The right-hand 
panel contrasts the results for p = 0.75 (brown stars) and p = 0.25 (blue squares). 
Whenp = 0.75, the threshold strength Z c -|_(£) (brown stars, dashed line) increases 
with £ while Z c _(£) (brown stars, solid line) remains approximately constant. 
This surprising result implies that the broken VT symmetry at Z > 4.5, !; = 
is restored by adding localized-gain £ ~ 5 to the constant-gain region (dashed 
arrow) . 

Now we present the phase diagram in the (Z, £) plane for < £, < £ c (p) as a 
function of the localized-gain location p. The left-hand panel in Fig. [2] shows the 
phase boundaries Z c+ {^) (dashed lines) and Z c ^{£) (solid lines) for p — 0.5 (black 
circles), p — 0.25 (blue squares), and p — 0.20 (red diamonds). When £ = 0, the phase 
boundaries are symmetric about the vertical axis, Z c+ (£) = — Z c _(£) ~ 4.48 |15j . 
When £ > 0, the phase boundaries become asymmetrical about Z = 0. In general, 
as the gain strength £ > for the localized potential increases, the corresponding 
thresholds for the constant potential Z c decrease. This behavior is not unexpected. 
The 'PT-symmetry is unbroken when the "effective potential strength" , naively 
proportional to |£ + Z\, is smaller than a threshold value. Therefore, as £ increases, 
maintaining the same "effective strength" requires lowering the threshold values 
Z c ± (£) including driving them negative. 

As we will show now, this argument is not generally applicable. The right-hand 
panel in Fig.[2]shows the PT-phase boundaries for p = 0.75 (brown stars), along with 

f Note that for certain rational values of p including p = 0.5, the eigenvalue _E n +i(£) is independent 
of ^, and energies (S„,E n +2) become degenerate leading to "PT-symmetry breaking [18] . 
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results for p — 0.25 (blue diamonds) for comparison. When p = 0.75, as the localized 
gain strength £ increases from zero, the corresponding gain threshold Z c+ (£) for the 
constant potential strength also increases whereas the loss threshold Z c _(£) remains 
essentially constant. It implies that a broken 'PT-symmetric phase that occurs at, 
say Z > 4.5 and £ = 0, will be restored by the addition of a localized-gain impurity 
£ ~ 5 at position p = 0.75 to the constant-gain region. Thus, contrary to standard 
expectations, we predict that the PT-symmetric phase is strengthened by adding two 
gain potentials. We emphasize that this re-entrant behavior of the PT-symmetric 
phase is generic and becomes more prominent as p — > 1. This surprising result shows 
that the two complex, "PT-symmetric potentials can lead to a co-operative effect that 
enhances the PT-symmetric threshold. 

4. Continuum vs. lattice results 

This continuum model has nonzero thresholds Z c or £ c (p) below which the VT- 
symmetry is unbroken. On the other hand, the lattice models of these two potentials 
have a fragile PT-symmetric phase, defined by vanishing threshold strengths, except 
when the localized gain or loss impurities are closest or the farthest [THJ HOI HI] ■ We 
now focus on the apparent discrepancy between these two results. Let us consider a 
lattice with N 3> 1 sites, nearest-neighbor spacing a/, and tunneling to, with open 
boundary conditions. In the continuum limit, a^ — > 0, N — > oo such that 2L — Na^ 
remains finite. In addition, when a^ — > 0, the tunneling to — > oo such that the product 
toa| = ti 2 /2m remains finite and gives the "effective mass" m of the particle hopping 
on this lattice. 

For the localized gain and loss potentials the 5-function strength £ = a^7 where 
±17 denotes the on-site impurity potential in the lattice version. Thus, the continuum 
and lattice thresholds for fractional position p are related by £(p)/(h 2 /2mL) = 
Nj c (p)/to- For almost all values of p, the lattice phase is fragile, j c (p)/to oc N^ 1 — > 
as N — > oo [21]. However, due to its precise TV -1 scaling, the continuum threshold 
£ c (p), measured in units of continuum scale, remains nonzero [15] . On the other hand, 
finite threshold strength 7 c /to ~ 1 for closest or farthest impurities implies that the 
corresponding continuum threshold diverges as £ c /(fi 2 /2mL) ~ N. This is precisely 
the result shown in the right-hand panel in Fig. [lj closest or farthest impurities on a 
lattice translate into continuum potentials located at a distance S oc 2/N away from 
the origin or the boundary respectively. 

For the constant gain and loss potential Vz, a corresponding argument shows that 
Z c / (ti 2 /2mL 2 ) = N 2 T c /to where the lattice- version of Vz(x) is given by the on-site 
potential Vj — irsign(j — N/2) for site index j = 1,...,N, and r c is the critical 
impurity strength for the lattice Hamiltonian. Since the continuum threshold is finite 
Z c / (h 2 /2mL 2 ) ~ 4.48, we expect that the corresponding PT-symmetric phase on the 
lattice will be fragile with scaling behavior T c /to oc 7V~ 2 . Figure [3] shows the critical 
threshold T c (N)/to for lattice sizes N = 20 to N = 2500; note the logarithmic scale 
on both axes. The dotted horizontal and vertical lines span one and two decades 
respectively, and thus the numerically obtained results (blue circles) show excellent 
agreement with the prediction that follows from continuum-limit considerations. 
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Figure 3. PT-symmetry breaking threshold r c (A r )/to as a function of the lattice 
size N for a tight-binding lattice with tunneling to and complex, on-site potential 
Vj = irsign(j — N/2); note the logarithmic scale on both axes. Robust phase of 
the corresponding continuum problem I15I implies that the corresponding VT- 
phase on the lattice is fragile with the lattice threshold r c /to oc N~ 2 ; the results 
for lattices with N = 20 to N = 2500 (blue circles) show that it is indeed the 
case. 



5. Conclusions 

In this paper, we have investigated the competition between two complex, PT- 
symmetric potentials on the PT-symmetric phase diagram. We have shown that a 
(simple) system with such potentials supports an unusual phase diagram and the 
restoration of PT-symmetric phase. At present, the most promising or current 
candidates for experimental realization of a PT-system can be modeled by tight- 
binding lattices and are explored via time-evolution of an initial state \<p) that is 
localized to a few, if not single, lattice sites [SHU]- Such a localized state contains 
energy components across the entire lattice bandwidth 4i , and therefore can only 
probe the fragile PT-symmetric phase of the lattice. 

We have shown that the robust or fragile nature of the PT'-symmetric phase 
is essentially determined by the relevant, continuum or lattice, energy scales. Thus 
PT-symmetry breaking effects in the continuum problem, including re-emergence of 
the PT-symmetric phase with increasing potential strength, can be probed through 
the dynamics of a broad and smooth initial state |$) with size 1 <C M -C N that 
contains only the energy components near the bottom of the energy band. An 
experimental investigation of such a system will deepen our understanding of PT- 
symmetry breaking in continuum and lattice models. 
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